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Riassunto. Viene fornita una rassegna di modelli matematici, equazioni differenziali 
nonlineari di tipo ellittico-parabolico, nei processi di elaborazione di immagini e di sequenze 
di immagini. 

Vengono inoltre presentati risultati recenti e problemi ancora aperti che riguardano l’e- 
sistenza, la regolarità e la stabilità delle soluzioni di alcune equazioni modello. 


Summary. We present somo mathematical models, nonlinear differential equations of 
elliptic-parabolic type, arising in immage processing. 

We also present some recently results and open problems about the existence, the regu- 
larity and the stability of solutions of such equations. 
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1. Preliminari. 


Un’ immagine in bianco e nero si può rappresentare mediante una funzione u 
definita in RN (con N = 2 oppure N = 3) dove u(x) rappresenta il livello di grigio 
nel punto x. 

Il primo passo in un processo di elaborazione di immagini è quello di ridurre il 
rumore presente nell’immagine per ottenerne una versione modificata che contenga 
informazioni significative dell’ originale. L'immagine data uo è quindi trasformata 
in nuove immagini e il principale parametro di questa trasformazione è la ”scala” 
la quale fornisce il grado di regolarizzazione o in un certo senso, la misura dell’ 
intorno del punto x considerato per stimare il livello di grigio in £. 

L’analisi multiscala è quindi una famiglia di trasformazioni (T):>o che applicata 
alla figura originale up fornisce una sequenza di immagini semplificate 


u(r,t) =(Truo)(2), 


la funzione u(x,t) è detta analisi dell'immagine alla scala t. 

L. Alvarez, F.Guichard, P.L. Lions e J.M. Morel in [1] hanno dimostrato che 
sotto ragionevoli ipotesi su (T:):>o (principi formali o assiomi) la funzione u(2,t) è 
soluzione di un’equazione parabolica con dato iniziale l'immagine originale. 


Il classico modello multiscala formalizzato da Witkin nel 1983 in [20] è costituito 
dall’equazione del calore ed è l’unico modello in cui (T:):>o è lineare ed invariante 
per isometrie. Questo modello nasce dall’ idea di realizzare un’analisi multiscala 
che regolarizzi la figura iniziale uo mediante una convoluzione con una funzione 
G gaussiana G(x,t) = ct-!/?exp(—|z|?/(4t)). Il principale svantaggio di questa 
regolarizzazione è che non solo riduce il rumore ma sfuma anche i contorni della 
figura rendendo così difficile la sua identificazione. Più precisamente, la linearità 
dell’operatore produce una perdita di edges cioè dei profili delle figure presenti 
nell’immagine, ossia l’insieme dei punti in cui il gradiente spaziale dell’immagine 
regolarizzata ha un massimo. 

Un importante miglioramento del modello fu introdotto nel 1987 da.J. Malik e P. 
Perona in [15]. Essi proposero di sostituire l'equazione del calore con un'equazione 
non lineare dei mezzi porosi del tipo 


a) | duu= div(g(1Dul)Du), 


dove g è una funzione di classe C©®°(R) non negativa, non crescente, tale g(0) = 
1 e g tende a zero all’infinito. In questo modo il processo di regolarizzazione è 
condizionato, nel senso che, se il gradiente spaziale di u è grande allora la diffusione 
è piccola e quindi è possibile una precisa localizzazione degli edges. Se il gradiente 
in un punto è piccolo la diffusione tende a regolarizzare maggiormente vicino a quel 
punto. 

I risultati sperimentali furono incoraggianti (gli edges infatti rimangono stabili 
lungo la scala) tuttavia l’equazione presenta serie difficoltà teoriche. Infatti, possi- 
amo scrivere 


div(9(|Du|)Du) = g(|Du|)Au + g'(|Du|)|Du|7!((D°u)Du, Du) 
= g(|1Du|)ucc + H'(|Du|)unn 
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dove H(s) = sg(s), 
EE (D°u)Dut Du* 
€ Du [Dut] 


è il termine di diffusione nella direzione ortogonale al gradiente e 


dis (D°u)Du Du 
"\ Du] |Dul 


è il termine di diffusione nella direzione del gradiente. 

Poichè g è nonnegativa il termine g(|Du|)uc produce una diffusione nella di- 
rezione degli edges. 

Se la funzione H” (il coefficiente di diffusione nella direzione ortogonale agli 
edges) è non negativa allora l'equazione è di tipo parabolica e vale il principio del 
confronto. Se invece H'(s) è negativa per qualche s allora l'equazione (1) cor- 
risponde ad una equazione del calore retrograda e quindi non ci si può aspettare 
unicità della soluzione. (Si veda [12] per un semplice esempio). 


Nel 1992, L. Alvarez, P.L Lions e J. M. Morel in [2] hanno proposto un modello 
basato sull’equazione del moto per curvatura media 


(2) Ou = |Du|div(Du/|Dul). 


Questa equazione non produce soltanto una diffusione di u nella direzione orto- 
gonale al suo gradiente. Infatti, 


du = |Du|div(Du/|Du|) = Au- ((D°u)Du, Du)/|Du|? 


e il secondo termine corrisponde ad una inibizione della diffusione nella direzione 
del gradiente. (La funzione u viene regolarizzata maggiormente nei punti lontani 
dagli edges, con una regolarizzazione minima degli edges stessi). In questo modello 
gli edges non sono altro che i bordi degli insiemi di livello di u. 

L'equazione (2) ha una sua interpretazione geometrica: almeno formalmente gli 
insiemi di livello della soluzione u(-,t) evolvono nella direzione normale con una 
velocità proporzionale alla loro curvatura media, almeno nella regione in cui u è 
regolare e il gradiente di u è diverso da zero. Una teoria delle soluzioni deboli 
basata sulla nozione di soluzione viscosa (cfr. [8]) fu proposta da Y.G.Chen, Y. 
Giga e S.Goto in [7], L.C.Evans e J. Spruck in [9], e da molti altri autori. 

Successivamente, F. Cattè, P.L. Lions, J. M. Morel e T. Coll in [6] hanno con- 
siderato il seguente modello, che è una sintesi dell’idea di Perona e Malik, ma che 
è consistente dal punto di vista formale 


du = div(9(|DGx*u|)Du) in Qx]0, T[, 
(8) u(7,0) =uo(2) on Q, 
dove 

Go(x) = co? exp(—|z|? /40) 


e la funzione g è come nel modello di Perona-Malik. 

Poichè la funzione G(2,t) = G:(x) è la soluzione dell’equazione del calore, il 
termine (DG; + u)(r,t) è il gradiente, al tempo 0, della soluzione del calore con 
dato iniziale u(z,t). 
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In [6], viene provato che il problema (3), con dato iniziale uo € L?(9) e dato al 
bordo di tipo Neumann, ha un'unica soluzione di classe C9(Qx]0, T[). 

(Il dato al bordo di tipo Neumann è quello più naturale nel processo di elabo- 
razione di immagini, infatti, corrisponde ad una riflessione della figura attraverso 
il bordo e ha il vantaggio di non imporre nessun valore al bordo e quindi di non 
creare un edge su di esso). 


Infine, L. Alvarez, P.L. Lions e J.M. Morel hanno proposto un modello basato 
sull’equazione del moto per curvatura media e sul metodo di stima precedente 


du = g(|DG, + u|)|Du|div(Du/|Du|) in Qx]0,T[, 
(4) u(z,0) = uo(z) on L. 

In (4) il termine g(|DGg + u|)| controlla la velocità di diffusione: se il gradiente di u 
ha media piccola in un intorno di x allora x è considerato un punto interno di una 
regione regolare dell’immagine e quindi la diffusione è forte. Altrimenti, il punto x 
viene considerato un punto di edge e la velocità di diffusione viene diminuita. 

In [2] viene provato un teorema di esistenza e unicità (in senso viscoso) per 
soluzioni continue in RN x [0,+00[, appartenenti ad L°(0,T;W!°(RN)), per 
ogni T > 0, (W!*(R”) è lo spazio delle funzioni lipschitziane e limitate), qualora 
Uo € W1°(RN). 

Inoltre se u e v sono soluzioni del problema (4), con dati iniziali rispettivamente 
uo € vo, allora esiste una costante positiva c tale che 


sup ||u(-,t) — v(-,t)]lz>(aN) £ c]luo — vollz®(r»)- 
1€[0,T] 


2. Assiomi ed equazioni fondamentali nell’elaborazione di immagini 


Come anticipato nell’introduzione, L. Alvarez, F. Guichard, P.L. Lions e J.M. 
Morel in [1] hanno dimostrato che le analisi multiscala (T;):>o sono governate in 
modo naturale da equazioni di tipo parabolico con dato iniziale l’immagine originale 
e queste equazioni possono essere dedotte da assiomi naturali sull’analisi multiscala. 

Il modello più semplice è quello della piramide visuale il cui dato iniziale è la 
figura originale uo, l'output è l’immagine multiscala u(7,t) e t è un parametro 
che può essere identificato con la misura dell’ intorno di x che è stato considerato 
per determinare il valore di T;uo(r). La successioni di immagini è così ottenuta 
mediante una sequenza di’ filtri (T.)>o sulla quale vengono fatte delle richieste 
formali che si possono distinguere in assiomi architetturali (causalità, regolarità), 
assiomi di stabilità (principio del confronto) e assiomi morfologici (invarianza per 
traslazioni, invarianza di scala). 


Assiomi architetturali. 
[Causalità] (L'immagine alla scala t può essere dedotta dalle immagini alle scale 
precedenti). 
Esiste una famiglia di operatori di traslazione (T{,3)1,s>o tali che To = 1 e 


Ti+3=T1+3,30Ts Vs,t,h>O0. 


[Località] (Il valore di T,f per t piccolo, in ogni punto 7, è determinato dal 
comportamento di f vicino ad x). 
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i (T..sf-—T,sg)(r)=o(t-s) pert-s-+0t, 
per tutte le funzioni f,g € Cf° tali che Df,Dg #0e D°f=D°gin RN, Vla|> 0. 
[Regolarità] 
IlTt+a,3(f +h9) — Ti+s,sf — hglloo <est Vs,t,h € [0,1] 
e Vf,g € C}° tali che Df,Dg #0. Inoltre 
||Te+s,ef — Tsoflloo < esn(t) 


e Vf € C}° tali che Df #0. Dove n(t) + 0 se t tende a zero, uniformemente in 
s €]0, 1]. 
Principio del confronto. 


[Principio del confronto] (Se un’immagine è più chiara di un’altra di un’altra 
allora quest’ordine si conserva lungo l’analisi). 
T.(f)<T:(9) inR", Vt>O 


e Vf,g € Cg° tali che Df,Dg #0. 
Assiomi morfologici. 


[Invarianza di scala di grigio] (L’analisi non si altera se cambia il contrasto). 


Ti(h(f))=h(T:(f)) per tutte le funzioni è monotone crescenti 
e Vf € C5° tale che Df#0. 
[Invarianza per traslazione] 
Ti(ta(f)) = ra(T:(f)) VheRN,t20 vfeCP,Df #0, 
dove rrf(7) = f(cr +h). 
Teorema (Theorem 2 in [1]). Supponiamo che (Ti): soddisfi gli assiomi prece- 


denti. Allora esiste una funzione F = F(A,p,x), definita su SN x RN xRe 
continua per p# 0, tale che 


(i) situ) e f 


uniformemente su RN (escluso i punti in cui Df = 0) e Vf € Cp° tale che Df # 0; 
(ii) F(A,p,x) > F(B,p,r) VpeRN,t>0, VA, BE SNcon A> B. 


(F è un generatore infinitesimale del semigruppo T\). 
Inoltre, se uo è una funzione uniformemente continua e limitata, allora la fun- 
zione 


+ F(D°f,Df,t) pers-+0 evi>0 


u(2,t) = T;(uo)(2) 
è l’unica soluzione viscosa continua e limitata su RN x [0,+00[, uniformemente 
continua in x, uniformemente rispetto a t, del problema 


Ou = F(D°u, Du,t) .in RN x [0,+00], 
u(z,0)=uo(7) in RN. 


(5) 


Questo approccio all'analisi multiscala consente di caratterizzare assiomatica- 
mente la maggior parte dei modelli già presenti e di introdurne dei nuovi. 
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Teorema. Supponiamo che (T:): verifichi gli assiomi precedenti escluso l’invarian- 
za per scala di grigio. Supponiamo inoltre che T; sia lineare, invariante per isome- 
trie, cioè 


T,.s(Rf)= RTs(f) Vi,s>0,vVf € CH, Df #0 
dove (Rf)(x) = f(Rx), R è una trasformazione ortogonale di RN e 


Ti(f+c)=T:(f)+c per ogni costante c. 
Allora ux,t)= Ti(uo)(x)è soluzione dell’equazione del calore 


du=pAu in È x [0,+00], 
u(z,0) = vo(r) in RN. 


dove pu è una opportuna costante positiva. 


Osservazione. Il modello di diffusione anisotropo di Perona-Malik soddisfa gli 
assiomi architetturali e le elementari proprietà di invarianza (isometria e traslazioni 
nello spazio). L’invarianza per scala di grigio non è mai soddisfatta. Infatti, se lo 
fosse allora F(XA, Ap) = AF(A,p), VA > 0 e VA, p, e poichè 


F(A,p) = g'(lp|)Ip| (Ap, p) + g(Ipl) traccia(A) 


scegliendo la matrice simmetrica A con traccia diversa da zero e tale che Ap = 0, 
si ottiene g(A|p|) = g(lp|) per ogni A > 0 e quindi g è necessariamente una funzione 
costante. 

Anche l’equazione del calore non la soddisfa, infatti se u è tale che Qu = Aue 
u= h(v), allora h'(v)0.v =  h'(v)Av +h"(v)|Dv|? e quindi h"(v)|Dv|? = 0 e poichè 
v è arbitrario si ottiene h" = 0. Concludendo, l'equazione del calore è invariante 
solo per trasformazioni affini della scala di grigio. 


Teorema. Se (T;): è un'analisi multiscala che soddisfa l'assioma di invarianza per 
di scala di grigio e se N =2, allora la funzione F in (5) assume la forma seguente 


(6) F(D?u, Du,t) = |Du|G (curo). Taj Dul' 3) 


dove la funzione G è una funzione continua in R x (R? \ {0}) x R,, nondecrescente 
rispetto al suo primo argomento, e curu(u)(x,t) è la curvatura dell'insieme di livello 
di u passante per x 


Du 
[Dul' 


curv(u) = div — 


Osservazioni. 
(1) Si osservi che se G = 1 (oppure G 
degenere che si ottiene è 


—1) allora l’equazione parabolica 


ul 


du=|Du| (oppure dQu= —|Dul). 
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Questa equazione è anche il modello fondamentale della matematica morfologica, 
d’altra parte dilatazione ed erosione possono essere interpretate come un moto di 
curve di livello con velocità costante nella direzione del gradiente. 

(2) Se l’analisi multiscala è anche invariante per isometrie e per trasformazioni 
affini allora la funzione G dipende solo dalla curvatura ed un modello possibile è il 
seguente 


(7) du = |Du|(curv(u))!/3. 


Questa equazione è stata trovata indipendentemente da G. Sapiro e A. Tannen- 
baum, (si veda [17]), nell’ambito dello studio dell’evoluzione di curve nel piano, 
come equazione del calore intrinseca 


(8) 3 Ci = Css 

dove s è la lunghezza d’ arco affine. (Questa equazione descrive il moto delle curve di 

livello di una soluzione di (7): è in questo senso che i due modelli sono equivalenti). 
Infatti, sia C : 5! x [0,T[+ R? una famiglia di curve regolari , C = C(p,t) con p 

parametro della curva. Supponiamo che la famiglia evolva secondo l’equazione (8), 

dove s è la lunghezza affine definita implicitamente dalla relazione 


det [C,, Css} = 1. 


Poichè 


(vettore normale affine) e Cp = &N con k, N, rispettivamente la curvatura media 
e la normale alla curva, si ha 


1 = det[C,, Css] = det[Cy, Oy] (2) =#(P\ 
= de 8) Uss] = de prpp de la da . 
Allora, nei punti in cui la curvatura è diversa da zero, si ha 


dp _psyy Pr 


Css = Cppk 2/9 + Cig ds? 


e quindi l'equazione (8) diventa 
d°p 
= p3N 4 Pr 
Ci= k/°N + 13? T 


2 x . . 
La componente tangenziale di può essere trascurata, infatti essa produce solo 
un moto dei punti lungo la curva stessa e quindi non ha nessuna influenza sulla 
evoluzione della curva intera. 
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In [16] viene provato che l’evoluzione secondo (8) conserva la convessità ed la 
famiglia di curve evolve ad un punto (diventando prima un ellisse), più precisa- 
mente: se Co è una curva regolare allora esiste una famiglia C = C(s,t) di curve 
regolari tali che 

aC = Css in St x[0,T[, 
Cli=o = Co. 


Inoltre, esiste to € [0, T[ tale che per ogni t > to la curva C(-,t) è convessa. 
Guichard in [11] ha invece provato che se il dato iniziale uo è una funzione 
lipschitziana allora l’equazione (7) ha un'unica soluzione viscosa u(£, t) lipschitziana 
in x e hélderiana in t di esponente 3/4. 
Osserviamo infine che l’equazione (2) del moto per curvatura media corrisponde 
all’equazione intrinseca Ci = C4s con s ascissa curvilinea. 


3. L’analisi di immagini in movimento. 


Un filmato oppure una sequenza di immagini può essere formalizzato mediante 
una funzione limitata u = u(2,0), dove 2 € R? è la variabile spaziale, 0 il tempo e 
u= u(2,0) il livello di grigio del punto 2 al tempo @. 

Consideriamo un’analisi multiscala del tipo 


du= F(D°u,Du,t), u=u(t,2,0) 


dove la funzione F = F(A, p,t) è non decrescente rispetto al primo argomento ed è 
continua se p # 0. 
Gli assiomi morfologici nelle variabili spaziali impongono in particolare che 


F(AA,Ap,t) = AF(A,p,t), VA>O, VA E S°, vpeR',p#0, 


F(A,p,t)= F(QA4Qp,p,t), perogni A€S*,pe R°,p70,Qp =In- PE 


Si possono inoltre introdurre nuovi assiomi che tengono conto della particolare 
natura del film. 


[Invarianza gallileiana] (Un moto con velocità costante dell'intera figura non 
altera l’analisi). . 

Se By(x,y,9) = (r — v19,y — 020,6) rappresenta un moto uniforme con velocità 
v= (01,02) allora 


F(BIAB,,Bip,t)=F(A,p,t), WeR?. 


[Invarianza rispetto al tempo] (L'analisi è invariante se il film viene uniforme- 
mente accelerato oppure decelerato). 

VI € Ret >0 esiste #' > 0 tale che la funzione t' = #'(A,t) è differenziabile 
rispetto À e 


: ar 
F(Sx\ASj, Syp,t) = pil (4,23), 


dove Sx(x,y,0) = (2,4, A0). 
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L. Alvarez, F. Guichard, P.L. Lions e J.M. Morel in [1] hanno dimostrato che 


l’unico modello dell'analisi multiscala verificante gli assiomi precedenti è il seguente 
(9) du = |Vu] curv(u) 3" (sgn(curv(u))accel?(u))*, ge [0, 1], 
dove Vu è il gradiente spaziale, curv è la curvatura spaziale delle linee di livello, 

\ 


((D?u)Du*, Du*) 


pi Put = (-upnuno) 


GRU) — 


e accel è un operatore differenziale del secondo ordine, detto accelerazione appar- 
ente, l’unico termine dell'equazione che fa intervenire le derivate rispetto al tempo 
0, definita nel seguente modo 


2 
2 Y 1 
(10) accel(u) = (x = 2) [Vul 


dove 


n = ((D°u)Dut, Dut), 
f2= ((D?u)Dut, Du), l 
%3 = ((D°u)Du*, Du*) 

e Dut = (uzuo,uyuo, —(u? + ui)). 

Significato dell’accelerazione apparente. 

Supponiamo che una figura si muova con moto di traslazione e che la funzione 
v= v(2,4,0) = (vz(0),vy(0)) rappresenti la velocità del punto (x,y) al tempo @. 
Supponiamo che il valore del livello di grigio di un punto non cambi durante il suo 
moto, cioè 

u(2(0), y(9), 0) = cost, 
allora derivando rispetto a @ si ottiene 
(0, Vu) + ug = 0. 


Siano vi la velocità nella direzione del gradiente spaziale 


_— (vu, Vu) _ uo 
015 Vul © [Vul 


e v2 la velocità nella direzione ortogonale al gradiente spaziale 


_ (V.Dut) 


7" Wa]! 


V=(v,1). 


Allora si può dimostrare che 


v2 = 12/(n1|Vul) 
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d 
(11) accel(u) = FG =-(Dv1,V), 
(la variazione di vi lungo la traiettoria del punto considerato). 
Supponiamo infatti che la figura w = w(x,y) si muova con moto di traslazione e 
u= u(x,y,0) sia la figura w in movimento, allora 


u(z,y,0)=w (- - f vz(s)ds,y — Î so) , 


Senza perdere di generalità possiamo supporre che al tempo 0 = 0 gli assi (î,j) 
siano tali che i sia nella direzione del gradiente spaziale e (5, Vu) = |Vu|. Allora 
Dw = (|Vw|,0,0), v1(0) = vz(0), v2(0) = vy(0). Inoltre 


Du = (|Vw|,0,-vi|Vul), 
Du = (0, |Vw|,0), 
Du* = |Vw|(-v1|Vw|,0,—|Vwl), 
è una base ortogonale di R?. 
Scrivendo poi la matrice hessiana di u e calcolando le sue restrizioni alla base 
prima definita si ottiene una espressione esplicita di Y1, Y2, 93 e di v2 in funzione 
di questi. 


Infine utilizzando la definizione (10) si ottiene la (11). 
Si può anche dimostrare che il vettore 


1 mn 
Va —-(DuÈ - Du) Ù 
Vul? ( N 


che è uguale alla velocità reale nei moti di traslazione, è tale che 


((D?u)V,V) a((D°u)V V a 
(DUVVI yi ( ) iva] 


Vul VI "IVI 


accel(u) = 


Infine, si osservi che (E, ma) |Vu|7! è la curvatura della curva di livello di u 


nel piano generato da Vu e V. 


Osservazioni. L’equazione (9) si può scrivere in termini di curvatura della 
superficie si livello della u = u(2,9) e della curvatura media spaziale degli insiemi 
di livello passanti per 2, ad esempio per q = 1/4 l'equazione diventa 


(n° _ AA 
(12) du = Vuli = |Du| (G (u)) 


dove G(u)(2,9) è la curvatura gaussiana della superficie di livello di u passante per 
(2,0). 

L’equazione (12) è stata individuata, utilizzando metodi completamente diffe- 
renti, anche da P.J. Olver, G. Sapiro e A. Tannabaum in [14] come un'equazione 
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di evoluzione di superfici, la più semplice equazione che sia invariante rispetto alle 
trasformazioni affini. 

Recentemente, Andrews in [4] ha provato che le superfici convesse di classe C? 
che evolvono secondo l’equazione (12) convergono ad un punto ellissoideo. (Si veda 
anche l’articolo precedente di Tso, [19]). 


Caselles e Sbert, hanno studiato sotto vari aspetti l'equazione 
Ou = |Du|segno(k1)((ki ko)t)", u= u(7,t), te R? 
dove ki e k2 in x sono le curvature principali della superficie di livello di u passante 
per x. Per o = 1/4si ottiene la (12). In [5], viene dimostrato, (anche per equazioni 
in una forma più generale), che se uo è una funzione continua e costante per |r| 


grande, allora esiste un’unica soluzione viscosa , continua e costante per |r| grande. 


Recentemente, Mikula, Sarti e Sgallari hanno proposto per l’elaborazione di 
sequenze di immagini tridimensionali, il modello 


Qu = |Vul (sgn(curv(u)) accel(u))* div (9(IDG5 * u|)Du). 


In [18] vengono forniti esperimenti numerici incoraggianti. 
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